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Abstract 

In this paper, a two-dimensional cellular automaton(CA) associated 
with a two-dimensional Burgers equation is presented. The 2D Burgers 
equation is an integrable generalization of the well-known Burgers equa- 
tion, and is transformed into a 2D diffusion equation by the Cole-Hopf 
transformation. The CA is derived from the 2D Burgers equation by using 
the ultradiscrete method, which can transform dependent variables into 
discrete ones. Some exact solutions of the CA, such as shock wave solu- 
tions, are studied in detail. 



1 



1 Introduction 



Recently, cellular automaton(CA) is extensively studied as a method to make 
models of complex systems, such as traffic flow[|^], fluid dynamics^, immune 
systems^]. Since CA models are represented by binary procedures or discrete 
equations, exact analysis or explicit solutions are rarely obtained. Thus we often 
use numerical calculations or thermodynamical and statistical methodsQ^ to 
study global features. 

Recently, the ultradiscrete method is proposed as a new type of discretization 
technique which allows us to study direct relations between continuous system 
and CA0[0. As for soliton equations, we can derive associated CA from a 
continuous equation keeping their integrable properties by the method. Another 
example of applying this method has been given by two of the authors and a CA 
associated with the ID Burgers equation is derived as follows^]: 

f/*+i = f/* + miniUU L-Ul)~ min(t/*, L - Ul,). (1) 

The CA is called Burgers CA(BCA), and it has nice properties; it can be trans- 
formed into a linear CA of diffusion type, and is a multi- value generalization 
of the rule-184 CAP|. Moreover, BCA and its extensions can be used as high- 
way traffic models, and properties of models and observed data are compared in 
detail|lOlin|. 



In this paper, we generalize the previous results to two spatial dimension. 
First, we introduce 2D Burgers equation, which can be linearized by a dependent 
variable transformation. Then we discretize independent variables keeping its 
integrable properties. Finally, we use the ultradiscrete method to obtain two- 
dimensional BCA. We also study some exact solutions of the CA. 

2 Two-dimensional Burgers equation and its 
ultradiscretization 

2D Burgers equation we consider in this paper is 

Ut = U^X + Uyy + 2UUa: + 2VUy, (2) 
Vx = Uy, (3) 

or, eliminating the variable v, we can rewrite these into a single equation 

.,^... + .„ + 2««. + 2«, (4) 

This equation has a remarkable property that it can be linearized by the Cole- 
Hopf transformation 

u=^ (5) 
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into a 2D diffusion equation 

ft = fxx + fyy (6) 

Tlius tlie Burgers equations are integrable in a sense tliat tliey can be transformed 
into a linear equation. It is noted tliat from (j^) and we obtain 

(7) 

From d^) and it is apparent tliat the variables u and v can be treated sym- 
metrically, then we can derive the time evolution equation for v. From (j^) and 
(1^, we obtain 

Vt = + Vyy + 2VVy + 2UV . (S) 

Thus we can write 2D Burgers equation in another coupled form as 

= V2u + 2(u- V)u, (9) 



where u = {u,v). This form appears in |12[, in which it is called vector Burgers 
equation. 

To discretize the independent variables in (||) and (H), we utihze discrete 
analogues to and Discretizing both time and space variables in (||), we 

obtain a discrete diffusion equation 

= (1 - 45)4 + ^(flij + fli^ + 4+1 + (10) 

where 6 = At/h'^ and h = Ax = Ay. At, Ax and Ay are lattice intervals in time 
t, space X and y, respectively. Next we define discrete analogues to the Cole-Hopf 
transformation and (0) 



<• - 4 - c% (11) 

J ij J tj 



where c is a constant. Evolution equations for ujj and vjj are derived from (|T0|) 
and (|lT]) as 

1-45 I I 1 I ^mj I 1 



" 1-45 <• 1 <■ 1 ' ^ ' 



(13) 



If we take an appropriate continuous limit At, h ^ 0, we obtain and (|^) 
from ( p!2D and (|13]). Moreover, (|12]) and (|13]) can be transformed into the linear 
diffusion equation (|TUp through (|TT]). 



Next, we ultradiscretize (0) and (|T3]), that is, discretize dependent variables 
u and V. Let us introduce a transformation of variables and parameters as follows: 

u% = e^/^ 4 = e^^/^ (14) 
= e--/^ (15) 



c5 

= e^/^ (16) 

where e is a parameter. Substituting these transformations into (0) and ([131) , 
taking a limit e — > +0 and using the formula 

A B 

lim e log(exp(— ) + exp( — )) = max(y4, B) = — min(— A, — -B), (17) 
£^+0 e e 

we obtain 

Ull' = f4 + min(M,f7t,^,L-f7j,V^5_„L-V^5) 

- min(M, f/*,, L - f/^^, L - VU,). (18) 

^^1,-^5 = f^+i-f^- (19) 

These are two-dimensional generalization of BCA ([^). We call this system two- 
dimensional Burgers CA(2DBCA). In (18), we can eliminate the variable V by 
using (^) and obtain 

f/j+i = f4 + g*,-g^+y, (20) 

where flow Q\j is given by 

Q% = min(M, [/^^, L - f/*,, \/o*^._i + E(f/*, - f/*,„i), L - V^^ - E(f/^,+i - 1/^,)), 

fc=0 k=0 

(21) 

which is an ultradiscrete version of (^. In (^I]) , V"oj- represents boundary condi- 
tions for at z = 0. Since variables m and v can be treated symmetrically, we 
can derive an equation for V by using a similar procedure to the above. The 
result is 

\^' = \% + mmiM,Ul,^,L-U!^,V^_,,L-\^) 

- min(M, L - V^, L - ^A^J. (22) 

It should be noted that 2DBCA is also expressed by a set of equations (ITS) and 
(p2D , and it is equivalent to that of ([T8|) and (|l^) under an appropriate condition. 
This fact is directly shown as follows. Let us deflne hjj by 

hl^ ^ min(M, f/*,, L - U^,^, VU,-i. L - VU,), (23) 
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then (0) and (p2D are written as 



u!;' = + K-i, - (24) 

V^;' = \^*+/iU-^U+i- (25) 



By using these equation we have 



T/t+l _ T/t+1 _ Tft+l I TTt + 1 

_ T/* _ r/* _|_ r/* 



^*H-^*-f/^,+i + f4- (26) 



Thus 



^Vi. - - f^+i + = const. (27) 

holds for any time t. By choosing initial condition of V and V so that the above 
constant becomes zero everywhere, then (|19D is always satisfied. Note that if 
U is periodic in lattice i and V in j, then the total sums Yl,i V\j and Yl,j are 
conserved quantities, which is directly shown by using (|18D and (|2^) . 

Finally, we consider a relation between the ultradiscrete Burgers equation and 
ultradiscrete diffusion equation. Introducing a transformation 

4. = exp(4/e), (28) 

ultradiscrete Cole-Hopf transformations 



U% = it,y-4 + -, (29) 
Vi = -4 + ^ (30) 



are obtained from (11) and ([T^ ) under the limit e — > +0. Substituting (^) into 
(|10|), we obtain an ultra-discrete diffusion equation; 

Fl+' = max(4 + § - M, i^^i,, ^-i,, " "^ax(0, § - M). (31) 

This is also obtained as follows. If we substitute (E9f) and (|30f) into (p!3), we have 



^i+lj ~ ™a^(4 ~'~ 2 ^i+lj^ ^i-lj^ 



- max(4i^. + - - M, 42^, 4, 4^+1, 4i,_i). (32) 



2" 

Since this equation is in a recurrence form in i, we obtain ( pi]) if we choose the 
decoupling constant in the above equation as — max(0, L/2—M). This shows that 
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we can linearize 2DBCA by the ultradiscrete Cole-Hopf transformations and the 
CA can be considered to keep integrable property. (|T9D is automatically satisfied 



by substituting (^) and (^) into it, and this means that (19) represents the 
compatibility condition of differences in i and j direction of the variable F. Note 
that we consider a linearity of (|31| ) in a sense that superposition of solutions is 
realized by max operation like '+' in (^. 

3 Some solutions of 2DBCA 
3.1 Shock wave solution 

In this section, we discuss some solutions and examine their behavior in detail. 
First, we derive an exact solution which represents shock waves. In order to 
obtain the solution, we utilize the shock wave solution of the discrete Burgers 
equation (|l^) and ([T3|) . Let us assume f-j has the following form 

N 

fij = l + T. ^Mkzt + r,3 + + e.), (33) 

2 = 1 

where k^, r^, ui^ and = 1, 2, ■ ■ ■ , A^) are constants. ( |55| ) is an exact solution 
of ([l0|) under the condition 

uj, = \og{l - A5 + 5{e^^ +e-^^ +e''^ +6-'^)) (;z = 1, 2, ■ ■ ■ , A^), (34) 

which is obtained simply substituting ( ^3]) into (p!0|). (0) is the dispersion relation 
between the frequency {uJz) and the wavenumber [k^ and r^). From (pAf ) we have 

_ ^ 1 + Ef=i exp(A;,(^ + 1) + + uj^t + ^^^^ 
1 + X;f=i exp(A;^i + r J + u^t + 



and also have v\y We call this 'W-shock wave" solution of (|T2p and (|T^). From 
this solution we can obtain a shock wave solution of 2DBCA by the ultradiscrete 
method. First we assume 

Kz Rz ^2 ^ fnn\ 

kz = — , = — , ujz = — , 4^ = — , (36) 

e e e e 

and noticing (|1^) and (p!6D, we obtain from ( p5| ) 

U^j = ^ + maxiO, 9i{i + 1, j, t), 62(1 + I, J, t),---,ON{i + I J, t)) 

- max(0, 6*1(2, j,t), 6*2(2, j,t),---, 6'Ar(i,j,t)), (37) 

by taking the limit e — > +0. Here phases 6z{z = 1, ■ ■ ■ , A^), are given by 

e,{t,j, t) = Kzi + Rzj + n.t + S,. (38) 
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From (0), a dispersion relation becomes in the limit 

= max(|K,|, \R,\, ^ - M) - max(0, ^ - M). (39) 

( p?!) is an exact solution of (|T8|) and we also get a solution by (|1^), which 
represent an ultradiscrete shock wave. The max operation in (0) represents the 
superposition principle, as mentioned in the previous section. 

Fig.l shows an example of shock wave solution (^) in the case of = 2. 
Parameters are taken as Ki = —1, Ri = K2 = -R2 = 1, = S2 = 0, L = 2 and 
M = 2. We easily see that hmf/*_. = L/2 + K2 = 2, and lim Uf. = L/2 + Ki = 0, 

«— »oo i— >— 00 ■' 

and lim f/*. = L/2 = 1. 

3.2 Particle model and excitation behavior 

If we assume the parameter M is sufficiently larger than L/2, we can neglect 
terms that contain M in (|1^), (0) and (|31|), and simplify the analysis on them. 
In this case equations on U and V become 

^1 = ^+mm{Ul,^,L-^,Vi_„L-\^) 

- min(t/^, L - Ul,^, VU^_,, L - V^^, (40) 

-min(^_i^.+„ L - Ul^^,, L - I/A+i), (41) 

and equation on F becomes 

i^^ = max(i^i^, ily, 4+1, 4„i). (42) 

It is easily verified that if initial conditions are taken as < 4, < Uf^ < L and 
< V^j < -Z^, then these conditions hold for any time t. 

The set of equations (^Oj) and (^Tj) express a simple particle model. Before 
explaining the model, we shortly review a one-dimensional particle model defined 
by (0) described in 0]. In the model, there are one-dimensional infinite sites 
indexed by i. Every site can hold L particles at most. Particles at each site move 
to empty spaces in its neighboring right site synchronously at each time step. 
According to this rule, fiow from site iioi + 1 at time t becomes min(t/*, L — U^_^i) 
if [/* denotes the number of particles at site i and time t. Thus considering the 
number of particles coming into and escaping from site i, we obtain (|l|). 



Next we explain two-dimensional particle model expressed by (|40|) and (^if ). 
The model is constructed as follows: 

(a) There are two-dimensional infinite sites indexed by i and j. 

(b) There are two kinds of particles, shortly saying U -particles and V^-particles. 
^/-particles move only in positive i direction. V^-particles in positive j di- 
rection, as shown in Fig. 2. 
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(c) Every site can hold L particles at most of both kinds respectively. 

(d) fZ-particles at site («, j) can move to empty spaces in its neighboring right 
site {i + 1, j) per unit time. V^-particles at (i, j) can move to those in up 
site (i, j + 1). 

(e) The number of [/-particles from {i — 1, j) to j) from t to t + 1 must be 
equal to that of V-particles from (i, j — 1) to (i, j) at the same time step. 

(f) Under this restriction, U- and V^-particles move to occupy neighboring 
empty spaces. 

Let us define the number of U- and V^-particles at site (i, j) and time t by Uj^ and 
respectively. Then, according to the above rule (e) and (f), flow of [/-particles 
from (i — l,j) to [i,]) and that of F-particles from — 1) to (i, j) are both 
Tcmi{Ul_^p L — Ulp V^j_i, L — Vlj). Thus, evolution equations of the above model 
become (^) and (^1]). We call this rule "pairing rule" because the same number 
of U- and l^-particles come into the same site. 

We show simple examples of time evolution of the above model. We choose 
L = 2 for all examples below. As the first example, we set = except 
F°o = 1 at the origin. This means U^^ q = 2, = 0, V^^^^ = 2, V^^^ = and 
f/°- = = 1 at all other points. The time evolution is given in Fig. 3. An exact 
solution of (|^) is 

^ij = ,I^J^t<t Gi+if^-k,j+jo-h (43) 

\k-l\<t 
k + l^t mod 2 

where Gk,i is defined by 

Gk,i = min(max(0,2A; 1) - 2max(0,2A;) max(0,2A; - 1), 

max(0,2/ + 1) - 2 max(0, 2/) + max(0, 2/ - 1)), (44) 

and io, jo is real phase constants. Since io and jo are arbitrary real constants, 
the above solution satisfies (^) even if i-j lattice is shifted as i ^ i + i'^ and 
j ^ j + Jo- If W6 set io = jo = in (^3D, we obtain the same solution shown in 
Fig. 3 on integer lattice points. 

We can derive the above solution from a solution of difference equation (p!0|) 
using ultradiscretization. The derivation is shown in Appendix. The solution 
represents an expansion of excited wave of F as shown in Fig.3(a). The solutions 
for U and V are obtained by substituting (|4^) into (p9D and (0), respectively, or 
simply by moving [/-particles and ^/-particles according to the above rule. Fig.4 
is an example of expanding rings of U and V. Initial conditions are Fj^j = except 
Fqq = 1 and F°o = 1- The solution is apparently obtained by a superposition of 
(^Sj). We see that Y.i U-j and Y.j ^Ij s-^e constant for time in both cases. 
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4 Concluding discussions 



In this paper, two-dimensional CA associated with Burgers equation(2DBCA) is 
presented by using ultradiscrete method. 2DBCA has some remarkable properties 
that it can be linearized by means of the ultradiscrete Cole-Hopf transformation, 
and has exact solutions which shows A^-shock wave and excitation behaviors. 

We also show a particle model of the CA, which can be considered as a pairing 
movement of particles in i and j direction. We hope that this rule and its inte- 
grability contribute to make models on some multi-dimensional complex systems. 
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Appendix 



In this appendix, we show that (^) can be derived from a solution of the differ- 
ence equation ( p!OD through the ultradiscretization. 

First, since we consider the case that the parameter M is sufficiently large, 
we can set 5 = 1/4 by using (0). Then ( [TI]| ) becomes 

fit^ — '^ifi+lj + fi~lj + fij+l + fij-l)- (^^) 

We can interpret this equation as 2D random walk process. Therefore, we obtain 
an exact solution including 2D binomial distribution, 

fij= Yl c\,i9i+io-K3+jo-h (46) 
ifc+ii<t 

\k-l\<t 
k+l — t mod 2 



where 



9k,i 



4t ^{t-\k + l\)/2)\{(t+\k + l\)/2)\ [{t-\k-l\)/2)\{{t+\k-l\)/2)\ 
(1 + exp(2A;/£))2 



+ exp((2A; + 1)/£))(1 + exp((2A; - l)/e)) 
^ (l + exp(2//e))2 



;i + exp((2/ + 1)/£))(1 + exp((2/ - 1)/^)) J 
and Jo is real phase constants, e is a limiting parameter. We have 



lim elog^ffc,/ = Gk,i in 

£— >+0 

and c\. I is always positive finite for finite /c, / and t. Therefore, taking a limit 
lini^e log all c\ i terms disappear and only gk^i survive in the summation of 

(^Bj). Thus we obtain F/^ in ( P3D from fj^ in by the ultradiscretization. 
Figure Al (a) shows F/^ , and (b) and (c) show e log fly We draw all figures as 
continuous plot on i and j and set t = 2 and = Jo = 0. Parameter e is 0.2 
and 0.01 in (b) and (c), respectively. From these figures, we can see a difference 
solution converges to a ultradiscrete one as e — > +0. 
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Figure Captions 

Fig.l The snapshot of shock wave solution at t = 1 in the case of = 2. Pa- 
rameters are Ki = —l,Ri = ^^^2 = -R2 = 1, 2i = S2 = 0, L = 2 and 
M = 2. 

Fig. 2 A physical explanation of the (2+l)-dimensional BCA, called "pairing rule" . 
The same number of fZ-particles and \^-particles come into the same site. 
The solid arrow and broken arrow show the movement of U and V, respec- 
tively. 

Fig. 3 Excitation behavior of (a)F, (b)f/, (c)^ at t = 6. We choose L = 2, and 
initial conditions are F^j = except the origin Fqq = 1. The light gray, 
dark gray and black square represent 0, 1, 2, respectively. 

Fig. 4 Expanding wave of U and V ((a)F, (b)f/, {c)V at t = 6). Initial conditions 
are F^ j = except Fq q = 1 and F^q = 1. The light gray, dark gray and 
black square represent 0, 1, 2, respectively. 

Fig.Al Plot of solutions of ultradiscrete and difference diffusion equation, (a) F/^ 
in (H), (b) ehgf^j in (H) for e = 0.2 and (c) ehgf^j in ® for e = 0.01. 
In all figures, we set t = 2 and io = jo = 0. 
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Figure 1 by K. Nishinari, et al 
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Figure 2 by K. Nishinari, et al 
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Figure 3 by K. Nishinari, et al 
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Figure 4 by K. Nishinari, et al 
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Figure Al by K. Nishinari, et al 
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